Monte Carlo studies of the chiral and spin orderings of the three-dimensional 

Heisenberg spin glass 
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The nature of the ordering of the three-dimensional isotropic Heisenberg spin glass with nearest- 
neighbor random Gaussian coupling is studied by extensive Monte Carlo simulations. Several in- 
dependent physical quantities are measured both for the spin and for the chirality, including the 
correlation-length ratio, the Binder ratio, the glass order parameter, the overlap distribution func- 
tion and the non-self-averageness parameter. By controlling the effect of the correction-to-scaling, 
we have obtained a numerical evidence for the occurrence of successive chiral-glass and spin-glass 
transitions at nonzero temperatures, Tea > Tsg > 0. Hence, the spin and the chirality are decou- 
pled in the ordering of the model. The chiral-glass exponents are estimated to be vcg = 1-4 ± 0.2 
and tjcg = 0.6 ± 0.2, indicating that the chiral-glass transition lies in a universality class differ- 
ent from that of the Ising spin glass. The possibility that the spin and chiral sectors undergo a 
simultaneous Kosterlitz-Thouless-type transition is ruled out. The chiral-glass state turns out to 
be non-self-averaging, possibly accompanying a one-step-like peculiar replica-symmetry breaking. 
Implications to the chirality scenario of experimental spin-glass transitions are discussed. 

PACS numbers: 



I. INTRODUCTION 

Spin glasses (SGs) are the type of random magnets pos- 
sessing both ferromagnetic and antiferromagnetic cou- 
plings, and are characterized by frustration and random- 
ness. The ordering of SG has been studied quite ex- 
tensively as a typical example of "complex system" and 
continues to give an impact on surrounding areas 
Most of theoretical and numerical studies on the SG or- 
dering have been based on a simplified statistical model 
called the Edwards- Anderson (EA) model, in which the 
spins are put on each site of a regular lattice and interact 
via the random coupling taking both positive (ferromag- 
netic) and negative (antiferromagnetic) signs [2| • For the 
Ising EA model in three dimensions (3D), it is now well 
established that the model exhibits an equilibrium SG 
transition at a finite temperature [H, 0, 0, S S H> @- 
The critical exponents of the SG transition evaluated by 
Monte Carlo (MC) simulations were compared favorably 
with those determined experimentally for the Ising-likc 
SG compound FeMnTi0 3 [TO]. 

Meanwhile, many of real SG magnets, including the 
well-studied canonical SGs, i.e., dilute magnetic alloys 
such as CuMn, AuFe and AgMn, are the Heisenberg- 
like magnets possessing only weak magnetic anisotropy. 
Thus, an isotropic Heisenberg EA model, rather than the 
strongly anisotropic Ising EA model, is expected to be a 
more realistic model. Experimentally, the existence of 
a finite-temperature SG transition and of a thermody- 
namic SG state in real Heisenberg-like SG material has 
been established 

For the 3D isotropic Heisenberg EA model, earlier nu- 
merical studies suggested, in apparent contrast to the 
experimental observation, that the model exhibited only 
a zero-temperature transition [lj, [l2|, [ID, Qj], [l5[ . Mean- 
while, one of present author (H.K.) suggested that the 
model might exhibit a finite-temperature transition in 



its chiral sector 0. Chirality is a multispin variable 
representing the handedness of the noncollinear or non- 
coplanar structures induced by frustration. It has sub- 
sequently been suggested that, in the ordering of the 3D 
Heisenberg SG, the chirality was "decoupled" from the 
spin, the chiral-glass (CG) order taking place at a tem- 
perature hi ghe r than the SG order, i.e., Tea > Tsg 
Ull H HM H3. Based on such a spin-chirality de- 
coupling picture of the 3D isotropic Heisenberg SG, a 
chirality scenario of experimental SG transition was pro- 
posed [ill, H3, HH [Hj]: According to this scenario, the 
chirality is a hidden order parameter of SG transition. 
Real SG transition of weakly anisotropic SG magnets is 
then a "disguised" CG transition, where the chirality is 
mixed into the spin sector via a weak random magnetic 
anisotropy. For a recent review, the reader is referred to 
Ref.Jp] 

The chirality scenario is capable of explaining sev- 
eral long-standing puzzles concerning the experimental 
SG transition [1] in a natural way, such as the origin 
of the non-Ising critical exponents observed in canonical 
SGs [20, HH , the apparent absence of the Heisenberg-to- 
Isin g cr ossover in the measured nonlinear susceptibility 
[2fl |2l|. and the origin of the mean-field-like transition 
lines (the so-called AT and GT lines) often observed ex- 
perimentally in the SG phase diagram in magnetic fields 
pol . I2H l23l . l24l [25j . The chirality scenario remains to be 
an attractive hypothesis in consistently explaining var- 
ious experimental observations for canonical SGs, and 
hence, it is an important task to examine the validity 
of the basic assumption underlying this scenario, i.e., 
the occurrence of the spin-chirality decoupling in the 3D 
isotropic Heisenberg SG. 

In recent numerical studies of the 3D Heisenberg EA 
model, consensus now seems to appear that the 3D 
Heisenberg SG indeed exhibits a finite-temperature tran- 
sition m F7L m, m, m m m m m h, h, m 
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IM H, H, H [H, However, the nature of the 

transition, especially whether the model really exhibits 
the spin-chiralit y d ecou pling, i s still under hot debate 
M, H, SM [H HI, H]. The present situ- 
ation is not completely satisfactory. Mentioning some 
of the recent numerical works: Hukushima and Kawa- 
mura studied the model with the random ±J coupling 
of modest lattice sizes L < 20 (L being the linear di- 
mension) but with a rather small number of samples of 
N s = 32 (for their largest L) These authors then 

concluded that T CG is certainly finite, while T$ G is either 
zero or nonzero but less than T G g, *- e -, T$g < Tcg, sup- 
porting the spin-chirality decoupling scenario. By con- 
trast, Lee and Young claimed on the basis of their data 
of the correlation- length ratio £/L of the model with the 
Gaussian coupling that the spin and chirality order at a 
common temperature, thus no spin-chirality decoupling 
[32l [3r| . However, their data suffer from either small 
lattice sizes of only L < 12 [32| or small number of sam- 
ples of N s — 56 [36| . Recently, Campos et al simulated 
the same model to much larger lattices L — 32 with 
larger number of samples iV s = 1000, but no data below 
the transition temperature [35]. Campos et al claimed 
that the chiral and spin sectors undergo simultaneously 
a Kosterlitz-Thouless (KT) transition with massive log- 
arithmic corrections. This interpretation, however, was 
criticized in Ref. [37| . 

Under such circumstances, we perform here a large- 
scale MC simulation of the 3D Heisenberg SG in order to 
shed further light on the nature of its spin and chirality 
ordering. We exceed the previous simulations by simulat- 
ing the system as large as L = 32 to temperatures consid- 
erably lower than T g for large number of samples of order 
N s ~ 10 3 . Note that none of the previous simulations 
satisfied all these criteria simultaneously. More impor- 
tantly, we calculate several independent physical quan- 
tities including the correlation-length ratios, the Binder 
ratios, the glass order parameters, the overlap distribu- 
tion functions and the non-self-averageness parameters, 
trying to draw consistent picture from these independent 
quantities, whereas Refs. [H, [H, Ull concentrated almost 
exclusively on the correlation-length ratio. By control- 
ling the correction-to-scaling effect in our data analysis, 
we can locate the chiral and spin transition points as 
T CG = 0.143 ± 0.003 and T SG = 0.125^;^. We then 
conclude that the SG transition occurs at a nonzero tem- 
perature which is located about 10^15% below the CG 
transition temperature. Thus, the 3D Heisenberg SG ex- 
hibits the spin-chirality decoupling. 

We also examine the possibility suggested in the previ- 
ous works that the spin and chiral sectors undergo simul- 
taneously a KT transition [35l [3r| . From our new data 
of large sizes L < 32 covering the temperature range be- 
low T g , we conclude that such a possibility can now be 
ruled out. In order to corroborate this conclusion, we 
also calculate the correlation-length ratio and the Binder 
ratio for the 2D ferromagnetic XY model, the standard 
model exhibiting the KT transition, and compare the re- 



sults with the ones of the 3D Heisenberg SG. Both the 
correlation-length ratio and the Binder ratio exhibit quite 
different behaviors between the two models, demonstrat- 
ing again that the transition of the 3D Heisenberg SG is 
not of KT-type. 

Recently, Pixley and Young questioned the utility of 
the Binder ratio in studying the ordering of vector SG 
models which is characterized by many-component ten- 
sorial order parameter 38]. To examine the validity 
of this claim, we also simulate the ferromagnetic 3D 
O(10) Heisenberg model, a model with a large number 
of order-parameter components n = 10. By calculat- 
ing the correlation-length ratio and the Binder ratio of 
the model, and by comparing them with those of the 3D 
Heisenberg SG, we conclude that the Binder ratio of the 
3D vector SG carries useful information independent of 
the correlation-length ratio, and should not be regarded 
as behaving in a trivial manner. 

We analyze the critical properties associated with the 
CG transition. On the basis of a finite-size scaling anal- 
ysis taking account of the leading correction-to-scaling, 
we get estimates of the CG critical exponents as vcg 
1 .4 ± 0.2 and r\cG = 0.6 ± 0.2, the former being the 
CG correlation-length exponent and the latter the CG 
critical-point-decay exponent. In fact, these CG expo- 
nents are close to the exponent values reported in earlier 
works, and also turn out to be very close to the SG expo- 
nent values experimentally observed for canonical SGs. 
In fact, this coincidence gives a strong support to the 
chirality scenario. 

In order to further probe the nature of the CG ordered 
state, we also calculate the overlap distribution function 
and the non-self-averaging parameter (the so-called A pa- 
rameter) both for the chirality and the spin. The CG 
ordered state turns out to be non-self-averaging. It also 
appears to exhibit a peculiar type of replica-symmetry 
breaking (RSB) which closely resembles the so-called one- 
step RSB. The behavior of the Binder ratio is fully con- 
sistent with such a one-step-like RSB picture. A prelim- 
inary account of the simulation was reported in Ref. (39j . 

The present paper is organized as follows. In §2, we 
define our model and explain some of the details of our 
numerical method employed. Particular attention is paid 
to the issue of thcrmalization, i.e., how we check the equi- 
libration which is often crucial in obtaining reliable data. 
Various physical quantities calculated in our simulations 
are introduced in §3. Then, our MC results are presented 
in §4. Quantities like the specific heat, the local-chirality 
amplitude, the CG and SG correlation-length ratios, the 
CG and SG susceptibility, the CG and SG Binder ra- 
tios, the CG and SG overlap distribution functions, the 
CG and SG non-self-averaging parameters etc, are calcu- 
lated. §5 is devoted to a finite-size scaling analysis of the 
CG critical properties. By analyzing the CG correlation- 
length ratio and the CG order parameter with taking 
account of the leading correction-to-scaling, we estimate 
the CG critical exponents. The character of the CG or- 
dered state is also studied via the Binder ratio, the over- 
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iV T 
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Tmin 
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2000 


32 


1 x 10 5 


0.333 


0.111 


8 


2000 


32 


1 x 10 5 


0.333 


0.111 


12 


2000 


32 


1 x 10 5 


0.333 


0.111 


16 


1500 


32 


1 x 10 s 


0.222 


0.121 


24 


1000 


44 


1 x 10 5 


0.222 


0.133 


32 


800 


48 


3 x 10 5 


0.209 


0.133 



TABLE I: Various parameters of our Monte Carlo simulations. 
L is the system size, N a is the number of samples, Nmc is the 
total number of Monte Carlo steps per spin (our unit Monte 
Carlo step consists of 1 heat-bath sweep and L over-relaxation 
sweeps), T m ax and T m i n are the highest and the lowest tem- 
peratures used in the temperature-exchange run, and Nt is 
the total number of temperature points. Measurements of 
physical quantities are made over the latter half of the total 
Nmc Monte Carlo steps, while the former half is discarded 
for thermalization. 



lap distribution function and the non-self-averageness pa- 
rameter. Experimental implications are briefly discussed. 
Finally, §6 is devoted to summary and discussion. 



II. THE MODEL AND METHOD 

We study an isotropic classical Heisenberg model on a 
3D simple-cubic lattice defined by the Hamiltonian 

H = - ^2 J v Si ' S i > W 
<ij> 

where Si — (Sf,Sf,S?) is a three-component unit vec- 
tor at the i-th site, and the < ij > sum is taken over 
all nearest-neighbor pairs. The couplings J^- are ran- 
dom Gaussian variables with zero mean and standard 
deviation unity. The lattice contains N — L 3 sites with 
L = 6, 8, 12, 16, 24, 32, periodic boundary conditions be- 
ing applied in all directions. 



We perform an equilibrium MC simulation by using the 
single-spin-flip heat-bath method and the over-relaxation 
method, which are combined with the temperature- 
exchange technique [icj . It has been demonstrated that 
this method is very effective in reducing the slow dynam- 
ics of hard-relaxing systems like SGs [361 ] . 

The simple cubic lattice consists of two interpenetrat- 
ing sublattices. We perform the heat-bath sweep sequen- 
tially through the sites on one sublattice after another. 
After the heat-bath sweep, we repeat the over-relaxation 
sweeps M times sequentially through the sites on each 
sublattice [4l|] . A unit over-relaxation process consists of 
computing the local field hi — J2j JijSj by a given 
spin Si and reflecting the spin Si with respect to the local 



field hi at this site. 

Si ► Si = Si + 2 — ~y£L — ~hi, (2) 

where hi = \hi\. 

The combination of one heat-bath sweep and M over- 
relaxation sweeps constitutes our unit MC step. In our 
following calculation, the number M is taken as being 
equal to the system size L, i.e., we take M = L. 

After every MC step, we perform the temperature- 
exchange trial. The method effectively promotes the 
system to overcome the free energy barrier characteris- 
tic of the spin-glass ordered state. We prepare Nt spin 
configurations with the same interaction coupling, some- 
times called "replicas" , which are located at distinct tem- 
peratures distributed in the temperature range between 
T m i n and T max . The maximum temperature T max needs 
to be high enough so that the auto-correlation time by 
the single-spin-flip dynamics is short enough. Then, the 
temperature-exchange trial is made between the two spin 
configurations at a pair of neighboring temperatures. 

In Table. I, we show some of the details of our sim- 
ulation conditions, including the system size (linear di- 
mension) L, the number of independent samples (bond 
realizations) N s , the number of temperature points used 
in the temperature-exchange process Nt, the minimum 
and maximum temperatures T m i n and T max , and the to- 
tal number of Monte Carlo steps per spin (MCS) per- 
formed per replica. The measurement is made over the 
last half of the Nmc MCS, while the former half is dis- 
carded for thermalization. The initial spin configuration 
is taken to be random. 

Error bars are estimated via sample-to-sample fluctua- 
tions for linear quantities like the order parameters, and 
by the jackknife method for non- linear quantities like the 
Binder ratio and the correlation length ratio. 

One of the most crucial issues in any equilibrium sim- 
ulation of SGs is to make sure that the system is fully 
thermalized. In particular, when we use an extended en- 
semble method like the temperature-exchange method, 
this point is particularly important, since, if only a part 
of replicas is not equilibrated, then all others might be 
affected, and the entire system might not be in equilib- 
rium. Hence, to ensure that the system is fully equili- 
brated, we need some stringent criteria for equilibration. 
In the present simulation, we have imposed the following 
six conditions for the check of equilibration. 
1) All of the "replicas" move back and forth many times 
along the temperature axis during the temperature- 
exchange process (typically more than 10 times) between 
the maximum and minimum temperature points. A typ- 
ical cycling pattern of a "replica" along the temperature 
axis during the temperature-exchange process is shown 
in Fig. 1(a) for our largest lattice size L = 32. We also 
check that the relaxation due to the single-spin-flip (1 
heat-bath sweep plus M = L over-relaxation sweeps) is 
fast enough at T = T max . Both the chiral and spin au- 
tocorrelation times at T = T max turn out to be about 22 



4 



0.209 



0.145 




3-10° 



0.112 



3-10° 



MC 



FIG. 1: (Color online) An example of typical "cycling" be- 
havior of a replica in the temperature-exchange run. The 
lattice size is L = 32. In (a), a replica exhibits a frequent 
cycling between T min and T max , where the minimum and 
maximum temperatures are chosen as T = T m i n = 0.133 and 
T = Tmax = 0.209, which correspond to the values of our final 
choice for L = 32. In (b), a replica exhibits a "trapping" be- 
havior, with its move limited in a narrow temperature range 
over long MC steps. In (b), we set the minimum temperature 
being lower, T m i n — 0.112, while the maximum temperature is 
the same T max = 0.209. Note that it often happens that, even 
when some of replicas exhibit a trapping behavior as shown 
in (b), other replicas exhibit apparently nice cycling behavior 
as shown in (a). When even a part of replicas exhibits such 
a "trapping" behavior, the system cannot be regard as being 
equilibrated. 



MCS for L = 32, and less than 16 MCS for smaller lattice 
sizes. This guarantees that different parts of the phase 
space are sampled in each "cycle" of the temperature- 
exchange process. 

It sometimes happens for larger lattice size L and the 
choice of lower T m i n that the frequent cycling between 
T m i n and T max cannot be achieved in a part of replicas: 
Some of the replicas, often not all, are "trapped" in a 
restricted temperature range in the course of simulation. 
A typical example of such a "trapping" behavior is shown 
in Fig. 1(b) for the case of our largest size L = 32, where 
the lowest temperature T m i n is taken to be T m i n = 0.112 
considerably lower than our final choice T m i n — 0.133. 
Once such a "trapping" occurs for certain replicas, an 
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FIG. 2: (Color online) Monte Carlo steps dependence of [qi] 
and [q s ] + ZL[E] defined by Eqs. (4) and (5). Thermal average 
is performed over the latter half of the total Nmc Monte Carlo 
steps. In equilibrium, these two quantities should coincide: 
See Eq. (3). Indeed, the curves of [qi] and of [q s ] + jjj{E] 
approach a common value within error bars. The lattice size 
is L = 32, and the temperature is T = T m i n = 0.133. The 
sample average is taken for a subset of total samples (150 
samples) . 



extremely long time is needed to get out of it, and the 
system can hardly reach thermal equilibrium. A partic- 
ularly tricky point here is that, once the trapping occurs 
for certain replicas, it often takes an extremely long time 
to get out of it so that various physical quantities appear 
to converge to "fake" stable values. It should also be 
noted that, even when certain replicas exhibit a trapping 
behavior as shown in Fig. 1 (b) , other replicas continue to 
exhibit a nice cycling behavior as shown in Fig. 1(a). Yet, 
if the trapping behavior is observed for a part of repli- 
cas, the system cannot be regarded as equilibrated at any 
temperature between T m i n and T max , since the ergodic- 
ity is not satisfied as an extended ensemble. Hence, we 
pay full attention that such a trapping does not occur 
in any replica by monitoring all the replicas. Once the 
trapping as in Fig. 1(b) is observed in certain replica at 
the time scale of, say, 10 5 MCS, we simply abandon the 
corresponding temperature set and try the new temper- 
ature set for all samples, not just one particular sample 
where we encountered trapping. 

2) We check that the relation expected to hold for the 
model with Gaussian bond distribution in equilibrium 
r42|, i.e., whether the relation 



2T 

h] = [<is} + j N [E} 



(3) 



is satisfied in the simulation. Here, E/N is the energy 
per spin and [• • • ] represents an average over the bond 
disorder. The "link spin overlap" qi and the quantity q s 
are defined by 



Ql 



(4) 



5 




MC 



FIG. 3: (Color online) The Monte Carlo steps Nnc depen- 
dence of various physical quantities, £cg/L and £sg/L (a), 
qcg and gsc (b), and q^Q and q^Q (c). Thermal average is 
performed over the latter half of the total Nmc Monte Carlo 
steps. The lattice size is L = 32, and the temperature is 
T = T m i n = 0.133. The sample average is taken for a subset 
of total samples (150 samples). 



in case of many other criteria. 

3) We check that measured physical quantities converge 
to stable values. As an example, we show in Fig. 3 the 
MC time dependence of several physical quantities, in- 
cluding the CG and SG correlation-length ratios £cg/L 
and £sg/L, the CG and SG Binder ratios gcG and gsc 

(2) (2) 

and the CG and SG glass order parameters q C Q and q S( j 
(to be defined in the next section) for L = 32 and at 
T = T m in = 0.133. All these quantities converge to 
stable values after some period, indicating that the sys- 
tem has been equilibrated. We note that the relaxation 
is faster at higher temperatures including the critical 
regime around Tqg- When equilibration is not sufficient, 
the correlation length £ and the glass order parameter 
q( 2 ) tend to be smaller than the true values as naturally 
anticipated, while no such inequality seems to exist for 
the Binder ratio g. As can be seen from Fig.3(b), it some- 
times occurs that thermalization of the chirality-related 
quantity takes more time than that of the spin-related 
quantity. This might indicate that the chirality is harder 
to relax than the spin. Hence, one should carefully test 
the stability of the chirality-related quantities in partic- 
ular, not only of the spin-related quantities. 

4) We check that the expected symmetry of the overlap 
distribution function Pj(q) under the reversal operation 
q — > —q holds for each individual sample. Since the global 
flipping of the spins and of the chiralities is supposed to 
give a slow mode of the system, this gives a stringent test 
of equilibration. Again, in our simulations, the symmetry 
of both spin and chiral P.j(q) turns out to be excellent 
for all individual samples [43|]. Of course, this is again 
only a necessary condition, not a sufficient condition. 

5) The equality between the specific heat computed via 
the energy fluctuation and the one computed via the tem- 
perature difference of the energy, which is expected to 
hold in any equilibrium system, is checked. 

6) We compare our data of the correlation length with the 
recent data reported by other authors, in the temperature 
range where common data are available [H, HH . 

We have carefully checked that our data satisfy all the 
criteria l)-6) above. In this way, we believe that the 
system has been fully equilibrated in our simulations up 
to the largest lattice size L = 32 and down to the lowest 
temperature T m i n . 



III. PHYSICAL QUANTITIES 



qs^il/N^d^-S,) 2 ), (5) 

ij 

where Nb — (z/2) is the number of nearest-neighbor 
bonds (z = 6 the coordination number of the lattice), 
and (• • • ) represents a thermal average. As illustrated in 
Fig. 2 in the case of our largest lattice L — 32, our data 
well satisfy Eq. (3) within the error bar. Note that, al- 
though this criterion is quite useful, it is only a necessary 
condition of equilibration, not a sufficient condition, as 



In this section, we define various physical quantities 
measured in our simulations and discuss some of their 
basics and details. 

For the Heisenberg spin, the local chirality at the i-th 
site and in the /i-th direction Xip, may be defined for three 
neighboring Heisenberg spins by a scalar 

Xin = S i+ e fi ■ (Si x Si-gJ, (6) 

where e M (/i = x, y, z) denotes a unit vector along the 
/i-th axis. There are in total 3Af local chiral variables. 
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First, we define an "overlap" for the chirality. In 
addition to "replicas" associated with the temperature- 
exchange process, we also prepare at each temperature 
two independent systems 1 and 2 (also called "replicas" 
here) described by the same Hamiltonian (1) with the 
same interaction set. We simulate these two replicas 1 
and 2 in parallel with using different spin initial condi- 
tions and different sequences of random numbers. 

The fc-dependent chiral overlap, <7 x (fc), is defined as an 
overlap variable between the two replicas 1 and 2 as a 
scalar 



f 



JV 



(1) (2) ik-f, 



(7) 



i—1 [i=x,y,z 



where the upper suffixes (1) and (2) denote the two repli- 
cas of the system. 

The fc-dependent spin overlap, q a p(k), is defined by a 
tensor variable between the a and (3 components of the 
Heisenberg spin, 



<7aM 



N 



(a,/3 = x,y,z). (8) 



In term of the fc-dependent overlap, the CG and SG 
order parameters are defined by the second moment of 
the overlap at a wavevector fc = 0, 



4g = [(?s(0) 2 ) 



j 2) _ [(|g x (Q)l 2 )] 

ICG - —4 



<7s(fc) 2 = ^2 k*/3(£) 

ot,/3=x,y,z 

(2) 



(9) 



(10) 



The CG order parameter q CG has been normalized here 
by the mean-square amplitude of the local chirality, 

1 



x 2 



3N 



N 



(ii) 



which remains nonzero only when the spin has a non- 
coplanar structure locally. The local-chirality amplitude 
depends on the temperature and the lattice size only 
weakly as shown later in Fig. 5. 

The CG and SG susceptibilities are defined by 



Xcg = SNq^l , xsg 
Finite-size correlation lengths are defined by 



Nqfl. 



1 



/ [(g(Q) 2 )] 

2sin(fc m /2)\/ [( g (fc m )2)] 



1, 



(12) 



(13) 



for each case of the chirality and the spin, £cg & n d Csg, 
where fc m = (27r/L,0,0) with fc m = |fc m |, and the \i- 
direction in Eq.(6) is taken here being parallel with fc. 



The CG and the SG Binder ratios are defined by 

1 / IM0) 4 )] \ 



9CG 



gsc 



3 - 




[(9x(0) 2 )] 5 



, [fe(o) 4 )] 
K? s (o) 2 )] 5 



(14) 



(15) 



These quantities are normalized so that, in the thermo- 
dynamic limit, they vanish in the high-temperature phase 
and gives unity in the non-degenerate ordered state. In 
the present Gaussian coupling model, the ground state 
is expected to be non-degenerate so that both gcc an d 
gsc should be unity at T = 0. 

One can also define the non-self-averageness parame- 
ter, or the so-called A parameter (44|, for the chirality 
and for the spin by 



Acg = 



[(^(o) 2 ) 2 ]-^) 2 )] 2 



2\12 



A S G = 



[<<?x(0) 



i(q s (d) 2 ) 2 } ~ [<9s(0) 2 >] 5 



[<9s(0) 2 >F 



(16) 



(17) 



The A parameter becomes nonzero if the CG or SG 
susceptibility is non-self-averaging. It should be noted 
here that, even if q S( j vanishes (or xsg remains finite), 
Asg could become nonzero if xsg is not self-averaging. 

Sometimes, one also uses the so-called Guerra param- 
eter, or the G parameter (45|, which is defined by 



G 



CG 



(9x(0) 



2\2l 



(9x(0) 2 )] 2 



[M0) 4 >] - [<? x (0) 2 )] : 



Gsg 



K<7s(0) 



2\2l 



K<7s(0) 2 )] 2 



[( 9S (0) 4 )]-[(<Z S (0) 2 



(18) 



(19) 



Unlike the A parameter, the G parameter can take a 
nonzero value even when the ordered state is a trivial one 
without accompanying an RSB [H, |47[. In fact, the G 
parameter is not independent of the Binder ratio g and 
the A parameter, given by 



Gcg — -^A C gI (1 ^ gcc), 



Gsg = 2 A sg/ (1 - gsc)- 



(20) 



(21) 



Hence, it should be noticed that, even if there is no SG 

(2) 

order in the sense q SG = and gsc — 0, Gsg could 
take a nonzero value if xsg is non-self-averaging, i.e., 
Asg^O. 

The chiral-overlap distribution P(q x ) is defined by 



P(q x ) = [(S(q x -q x (0)))}- 



(22) 
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The spin-overlap distribution P(qdiag) is defined origi- 
nally in the tensor space with 3x3 = 9 components. To 
make this quantity more easily visible, one may define 
the diagonal spin-ove rlap , which is a trace of the original 
tensor overlap as pjj [48| 



configuration is locally noncoplanar sustaining a nontriv- 
ial scalar chirality. 



P{qdiag) = [{5(qd 



lag 



E <w(o))>] 



(23) 



fi—x,y,z 



In the high-temperature phase, both P(q x ) and P{qdiag) 
should approach the ^-function at q = in the ther- 
modynamic limit. In the low-temperature phase, P{q x ) 
should develop two symmetric delta-function peaks at the 
q^-values corresponding to the chiral EA order parame- 
ter ±q^Q, while P(qdiag) should develop two symmetric 
delta-function peaks at 1 /3 of the spin EA order param- 
eter ±qsc : See Refs.[48j for further details. 



IV. 



MONTE CARLO RESULTS 



In this section, we present our Monte Carlo results on 
the three-dimensional isotropic Heisenberg SG with the 
random Gaussian coupling. 

We first show in Fig. 4 the temperature dependence of 
the specific heat for various lattice sizes L. An arrow 
in the figure indicates the location of the CG transition 
temperature Tcg, which will be determined below. As 
can be seen from the figure, the specific heat depends 
on the temperature only weakly without any appreciable 
anomaly. 
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FIG. 4: (Color online) The temperature and size dependence 
of the specific heat per spin. An arrow indicates the location 
of the chiral-glass transition point. The inset represents a 
magnified view. In the zero-temperature limit, the specific 
heat is expected to tend to unity. 

In Fig. 5, we show the mean-square local chirality am- 
plitude x 2 as defined by Eq.(ll). Due to the local na- 
ture of this quantity, it exhibits only very weak size de- 
pendence. It also depends on the temperature weakly, 
and tends to a nonzero value in the T — > limit, 
X(T = 0) ~ 0.274, indicating that the ordered-state spin 
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FIG. 5: (Color online) The temperature and size depen- 
dence of the local-chirality amplitude. An arrow indicates 
the location of the chiral-glass transition point. In the zero- 
temperature limit, the local-chirality amplitude is extrapo- 
lated to x 2 = 0.075 ± 0.001 (or X = 0.274 ± 0.002): See the 
broken line in the figure. 

The temperature dependence of the CG and SG or- 
der parameters q GG and q^ G are shown in Fig. 6 (a) and 
Fig. 6(b), respectively. The CG order parameter increases 
more sharply than the SG order parameter, suggesting 
that the chirality exhibits a stronger ordering tendency 
than the spin. Meanwhile, more careful analysis of the 
size dependence is required in determining the transition 
point, which will be postponed later in this section. 

The temperature dependence of the CG and SG sus- 
ceptibilities is shown in Fig. 7. In contrast to the SG 
susceptibility xsg which is found to be an increasing 
function of the lattice size L at all temperature studied, 
the CG susceptibility xcg behaves in this way only at 
T/J IS 0.165, but exhibits an opposite size-dependence 
at T/J > 0.165: See the inset of Fig.7(a). Since the 
chiral susceptibility in the critical regime should be an 
increasing function of L, this observation suggests that 
the critical region associated with the CG order might be 
rather narrow. Similar size dependence of XSG was also 
observed in an earlier work [2y|, and also in the 3D XY 

SG m. 

The temperature dependence of the CG and SG 
correlation-length ratios, £,cg/L and ^sg/L, is shown in 
Figs. 8 and 9, an overall behavior in Fig. 8 and an enlarged 
figure in Fig. 9. As can be seen from the figures, while the 
chiral £cg /L curves cross at temperatures which are only 
weakly L-dependent, the spin £,sg/L curves cross at pro- 
gressively lower temperatures as L increases. 

The present £/L data are compared with the data by 
other authors as follows: Our data for £/L are in full 
agreement with those of Ref. [35| within statistical er- 
ror bars over the narrow and relatively high-temperature 
range covered by their data. The data of Ref. [36( for their 
largest L (on which their claim for "marginal" was based) 
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FIG. 6: (Color online) The temperature and size dependence 
of the chiral- glass order parameter (a), and of the spin-glass 
order parameter (b). An arrow indicates the location of the 
chiral-glass transition point. Inset of Fig. (a) is an enlarged 
view of the transition region. 



are lower than our present ones and those of Ref . [35J by 
about 5 to 6 of our a units; this may be a purely sta- 
tistical effect in view of the limited number of samples 
measured in Ref. [3r| . 

As an other indicator of the transition, we show in 
Fig. 10 the Binder ratios for the chirality (a), and for the 
spin (b). The chiral Binder ratio gcG exhibits a nega- 
tive dip which deepens with increasing L. The data of 
different L cross on the negative side of gec- These fea- 
tures indicate a finite-temperature transition in the chiral 
sector. 

In order to estimate the bulk CG and SG transition 
temperatures quantitatively, we plot in Fig. 11 the cross- 
ing temperatures of £cg I L and £sg I L and those of £sg 
for pairs of successive L values versus 1/L av , where L av 
is a mean of the two sizes. The L at) -dependence of 
the dip temperature of gcG is also shown in the figure. 
Since the data turn out to show an almost linear 1/L av - 
dependence, we tried in Ref. [39j a simple linear extrapo- 
lation of the crossing temperatures T cross (L) and the dip 
temperature Tdi P {L) to obtain Tea — 0.145 ± 0.005 and 
T SG = 0.120 ± 0.006. In the following, we try further 
elaborate analysis. 

Considerable shift of the crossing temperature with the 



FIG. 7: (Color online) The temperature and size dependence 
of the chiral-glass susceptibility (a), and of the spin-glass sus- 
ceptibility (b). An arrow indicates the location of the chiral- 
glass transition point. Inset of Fig. (a) is an enlarged view 
of the temperature region somewhat higher than the critical 
regime. 

system size L observed in Fig. 11 suggests the relative 
importance of the correction-to-scaling effect. Generally, 
one expects 

T cross (L)-T cross (oo)^L- e , 9 = w + ^ (24) 

where v is the correlation-length exponent and to is 
the leading correction-to-scaling exponent. (Incidentally, 
Ref. [391] quoted 9 = to, which was inappropriate in the 
standard notation. This does not affect the subsequent 
analysis of Ref. [39], though.) Here we perform a joint 
fit of T cross (L) of both the chiral correlation- length ra- 
tio and the chiral Binder ratio, ticc/L and gec, to 
the form Eq.(24), where the CG transition temperature 
Tcg — T cross (oo) and the exponent 9 are taken to be 
common between S^cg/L and gec- The optimal fit is 
achieved at Tcg — 0.143 and 9 = 0.93. To estimate the 
error bar, we show in Fig. 12 the associated x 2 -values of 
the fit as a function of the assumed Tcg and 9 values: In 
Fig. 12(a), the dependence on Tcg is shown with optimiz- 
ing 9 for each Tcg, while, in Fig. 12(b), the dependence 
on 9 is shown with optimizing Tcg for each 9. From 
these plots, we get estimates Tcg — 0.143 ± 0.003 and 
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FIG. 8: (Color online) The temperature and size dependence 
of the correlation- length ratio for the chirality (a), and for the 
spin (b). 



9 = 0.93 ±0.06, which turn out to be consistent with our 
previous estimates T GG = 0.145 ± 0.005 and ~ 1 
Our present estimate of T GG also agrees with the one 
suggested by Campos et al, T GG — 0.147 (35[. The fact 
that we have two independent data sets for T cross (L), 
one from £, GG /L and the other from gcc facilitates our 
estimate of the chiral-glass transition temperature. 

We have performed a similar x 2_ analysis based on 
Eq.(24) also for the crossing temperature of the spin- 
glass correlation-length ratio £sc/L. (For the spin, wc 
have only one kind of crossing temperature.) We then get 
T SG = 0.125t^?2 and = 1.2 ±0.35. The obtained T SG 
value is consistent within the errors with our previous 
estimate T SG = 0.120 ± 0.006 [39|. 

Then, our estimates of T CG = 0.143 ± 0.003 and 
T SG = 0.125j£^ suggests that T SG is lower than T CG 
by about 10^15%, indicating the occurrence of the spin- 
chirality decoupling. If we force the chiral crossing- 
points data to obey T GG = 0.125, the associated x 2 - 
value is greater than the optimal value obtained with 
T GG — 0.143 by 25.5, which is significantly greater than 
the standard error-bar criterion, unity. Likewise, if wc 
force the spin crossing-points data to obey Ts G = 0.143, 
the associated x 2_v alue is greater than the optimal value 
obtained with T$ G = 0.125 by 19.3, which is again sig- 
nificantly greater than the standard error-bar criterion, 



FIG. 9: (Color online) Magnified view of the temperature 
and size dependence of the correlation-length ratio around 
the transition region for the chirality (a), and for the spin (b). 
An arrow indicates the location of the chiral-glass transition 
point. 

unity. Hence, a simultaneous spin and chiral transition 
is highly unlikely from our present data. 

In Fig. 13, we show the ratio of the CG and SG cor- 
relation lengths £cg / Lsg ■ The ratio curves of different 
L intersect. More precisely, for smaller sizes of L IS 12, 
this ratio tends to be almost size-independent at lower 
temperatures indicating that the chiral and spin correla- 
tion lengths behave quite similarly [32^ , By contrast, for 
larger sizes of L C 16, the ratio curve splays out in the 
lower temperature regime. This change of behavior of the 
ratio £cg / £sg is quite consistent with the expected size- 
crossover from the trivial coupling behavior for smaller 
L < 12 to the decoupling behavior for larger L > 16 [l9l ]. 
The crossing point of the ratio curves for our largest L 
comes around T ~ 0.154, which seems consistent with 
our estimate above of T GG ~ 0.143. Meanwhile, the ratio 
itself is still less than unity even at the lowest tempera- 
ture studied. Still larger lattice size is required to reach 
the region where this ratio exceeds unity. 

Recently, Campos et al claimed on the basis of their 
data of the correlation- length that the chiral and spin sec- 
tors undergo simultaneously a KT transition [35]. This 
interpretation, however, was criticized in Ref. [37j- More 
recently, Lee and Young also suggested on the basis of 
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FIG. 10: (Color online) The temperature and size dependence 
of the Binder ratio for the chirality (a), and for the spin (b). 
An arrow indicates the location of the chiral-glass transition 
point. 



their data of the correlation-length ratios that the sys- 
tem exhibits a "marginal" behavior [361 ] . In view of such 
recent claims on the model, we further examine here the 
possibility of a KT-type phase transition. 

First, we note that, with decreasing the temperature, 
our data of the CG and SG correlation-length ratios of 
various L do not merge as expected for the KT transition, 
but intersect, with crossing points shifting to lower tem- 
perature for large L. For temperatures below the crossing 
points, the £cg/L and £sg/L curves fan out, rather than 
becoming L- independent. In order to make a more direct 
comparison with the behavior of the KT transition, wc 
calculate the correlation-length ratio £/L for the ferro- 
magnetic 2D XY model, a standard model displaying 
the KT transition. The result is given in Fig. 21 of Ap- 
pendix A. As can be seen from the figure, £/L curves 
of the 2D XY ferromagnet do not cross at any temper- 
ature but merge for larger L, becoming asymptotically 
L-independent at temperatures lower than the KT tran- 
sition temperature. Hence, our present data of cither 
the CG or SG correlation-length ratio shown in Fig. 9 are 
radically different from the one of a typical KT tran- 
sition shown in Fig. 21: Our data of $ /L curves of the 
3D Heisenberg SG are not "merging" 35] nor "marginal" 
[3r3 |. but splay out. 
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FIG. 11: (Color online) The (inverse) size dependence of the 
crossing temperatures of £cg/L and £sg/L, the dip temper- 
ature Tdip and the crossing temperature T cr0 ss of gcG- Lines 
represent the fitting curves of the data based on Eq.(24). The 
spin-glass and chiral-glass transition temperatures are extrap- 
olated to T CG = 0.143 ± 0.003 and T SG = 0.125±g;g?|. The 
inset exhibits a wider range. 

In Ref . [35l| . Campos et al performed a KT-type scal- 
ing with massive logarithmic corrections for the SG 
correlation-length ratio £sg/L, and reported that the 
data exhibited a good scaling. We also tried an exactly 
same scaling plot with the same logarithmic correction 
term as performed by Campos et al (35| . and the result 
is shown in Fig. 13. Note that our present data include the 
low temperature range below Tqg which was not covered 
by the data by Campos et al. (The temperature range 
where Campos et al reported their scaling plot is indi- 
cated by the dashed-line box in our Fig. 13.) As is evident 
from Fig. 13, the KT scaling turns out to be poor, even 
with invoking a massive logarithmic correction. We thus 
conclude that the possibility that the spin and chiral sec- 
tors undergo a simultaneous KT-type transition can be 
ruled out from our present data of the correlation-length 
ratio. 

The peculiar form of gec with a negative dip shown 
in Fig. 10(a) is consistent with the occurrence of a one- 
step-like replica-symmetry breaking (RSB) as suggested 
by Hukushima and Kawamura [la , Il9| . This interpre- 
tation is corroborated by the form of the the calculated 
chiral-overlap distribution below Tcg to be shown later 
in Fig. 16 (a), which exhibits a prominent central peak at 
q x =0. 

By contrast, the corresponding spin Binder ratio gsG 
shown in Fig. 10(b) does not exhibit a crossing nor a 
merging in the temperature range studied, suggesting 
that the SG transition temperature, if any, occurs be- 
low T ~ 0.13. Meanwhile, as the size L is increased, gsG 
develops more and more singular form at low tempera- 
ture, indicating that the associated overlap distribution 
significantly changes its shape at low temperature. If one 
recalls the fact that gsc takes a value unity at T = 0, 
gsG is expected to develop a negative dip at a lower T 
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FIG. 12: (Color online) The x 2 -value of the fit of T cross (L) to 
the form Eq.(24) as a function of the assumed Tec-value (a), 
and of the assumed 8- value (b), where other fitting parameters 
are optimized. The best x 2 -value is obtained at Tec = 0.143 
and 6 = 0.93. The horizontal straight line corresponds to the 
X 2 -value equal to the minimum ^-vahie plus unity, which is 
usually used to estimate error bars. 



(of IS 0.13) accompanied by an upturn toward T = 0. 
This feature of gsG strongly suggests the occurrence of a 
SG transition at a nonzero temperature, T$g ~ 0.13. In 
order to locate T$g from gsc more directly, however, we 
need the large-lattice data at lower temperatures. 

In view of the suggestion of the KT-type transition in 
Refs. [35| and (36|, we further examine the possibility of 
the KT-type transition via the CG and SG Binder ra- 
tios. In order to make a direct comparison with a typical 
behavior of the KT transition, wc calculate the Binder 
ratio g for the ferromagnetic 2D XY model, and the re- 
sult is given in Fig. 22 of Appendix A. As can be seen 
from the figure, the g curves of the 2D XY ferromagnet 
of different L weakly cross at a temperature above the 
KT transition temperature Tkt, which gradually tend 
to Tkt in the L — > 00 limit. Hence, the present data 
of either the CG or SG Binder ratio shown in Fig. 10 are 
radically different from those of a typical KT transition: 
Our spin gsc curves do not cross, while our chiral gcG 
curves exhibit a negative dip. Again, our data of the 
Binder ratio of the 3D Heisenberg SG are inconsistent 
with the KT scenario. 
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FIG. 13: (Color online) The temperature and size depen- 
dence of the ratio of the CG correlation-length ratio to the 
SG correlation-length ratio for various sizes. An arrow indi- 
cates the location of the chiral-glass transition point. 
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FIG. 14: (Color online) Kosterlitz-Thouless-type scaling plot 
with a logarithmic correction as performed in Ref. 35] is ap- 
plied to our data of the spin correlation-length ratio. An ar- 
row indicates the location of the chiral-glass transition point. 
The window in the figure exhibits the data range presented 
in Fig. 3b of Ref. [35| . The scaling turns out to be poor. 



Pixley and Young recently criticized that the Binder 
ratio might not be an appropriate quantity in studying 
the ordering of vector SGs, arguing that the large num- 
ber of order-parameter components (n = 3 2 = 9 in the 
Heisenberg SG) might lead to a trivial Gaussian distribu- 
tion even below T g [38] . To check the validity of such an 
expectation, we calculate the Binder ratio g of a simple 
3D 0(n) ferromagnet with large number of n = 10 com- 
ponents, and the result is given in Fig. 23 of Appendix B. 
As can be seen from the figure, g of 3D O(10) Heisen- 
berg ferromagnet exhibits a clear crossing behavior at the 
transition temperature T c and splay out below T c , the be- 
havior characteristic of the standard long-range ordered 
phase, in spite of the large number of its order-parameter 
components. Very much similar behavior was also ob- 
served in the Binder ratio of the 3D 0(6) ferromagnet 
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[501 ] . Such a behavior of g is quite different from the one of 
the 3D Heisenberg SG we observed in Fig. 10(b). Hence, 
the result presents counter-examples to the criticism of 
Ref. [38| , demonstrating that the peculiar behavior of gsc 
observed for the 3D Heisenberg SG in Fig. 10(b) should 
be regarded as a manifestation of essential features of 
the SG ordering, not mere an artifact due to the large 
number of order-parameter components. 

In Fig. 15, we show the size dependence of the CG and 

(2) (2) 

SG order parameters q GG and q SG on a log-log plot for 
several temperatures. Straight lines are drawn by fit- 
ting the three data points of smaller sizes L = 6, 8 and 
12 at each temperature. As can be seen from Fig. 14(a), 

(2) 

q GG exhibits an almost linear behavior at a temperature 
T = 0.148, an upward curvature characteristic of a long- 
rage ordered state at lower T, and a downward curvature 
at higher T which should eventually tend to a linear be- 
havior with a slope equal to — d = —3 in the disordered 

(2) 

phase. Thus, the data of q GG are consistent with our 
conclusion from the analysis of £,cg/L and gec above 
that the CG transition occurs at T C g = 0.143 ± 0.003. 

(2) 

The SG order parameter q SG exhibits a significantly 
different behavior, i.e., it exhibits a downward curvature 
characteristic of a disordered state at T = 0.148 ~ Tcg, 
or even at T — 0.133 < Tcg- At our lowest temperature 
T/J= 0.121 where we could equilibrate only smaller lat- 
tices of L < 16, the data exhibit a near linear behavior up 
to L = 16, although it is not clear whether this linear be- 

(2) 

havior extends to larger L. Thus, our data of q SG (L) are 
consistent with our conclusion from the analysis of £sg / L 
above that a SG transition occurs at T$g — 0-125T^oi2> 

(2) 

whereas, from the present data of q SG only, we cannot 
rule out the possibility that T$g is significantly lower 
than this. Although reliable estimate of the correspond- 
ing SG exponents is difficult due the remaining uncer- 

(2) 

tainty in Tsc, °ur data of q SG in Fig. 14(b) enable us to 
conclude rise ~ —0.30, which definitely differs from the 
CG "t\cg value. 

In Figs. 16(a) and (b), we show the chiral-overlap distri- 
bution function (a) and the diagonal-spin-overlap distri- 
bution function (b) at a temperature T — 0.133 which lies 
below Tcg but above T$g- As one can see from Fig. 16(a), 
the chiral-overlap distribution function P{q x ) displays a 
central peak at q x = for L > 12 which grows with in- 
creasing L. In addition to the central peak, there exist 
small side peaks located at the g-values corresponding 
to the CG EA order parameter ±q§ G , though these side 
peaks are weak and look like "shoulders" at this temper- 
ature. At temperatures higher than Tcg, P(lx) exhibits 
only a single Gaussian peak at q x — without "shoul- 
ders" even for smaller lattices. The behavior of P(q x ) 
observed here is similar to the one reported before for 
the 3D Heisenberg SG with the Gaussian coupling [l~8| , 
the 3D Heisenberg SG with the binary coupling [19| and 
related Heisenberg SG models [H, |48[ . We note that the 
side peaks of P(q x ) were more clearly visible in Refs. [l8T] 
and [19( | . The form of the overlap distribution character- 
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FIG. 15: (Color online) The size dependence of the chiral- 
glass order parameter q^ a (a), and of the spin-glass order 
parameter q^Q (b). Straight lines in the figures are drawn by 
fitting the three data points of smaller sizes, L — 6, 8 and 12. 

ized by a central peak coexisting with side peaks is the 
one common to systems exhibiting the so-called one-step 
RSB. The observed feature of P(q x ) is also consistent 
with the existence of a negative dip in the CG Binder 
parameter gec and with the crossing of gec occurring 
on the negative side as discussed before. We note that 
a one-step feature was also suggested from the study of 
the fluctuation-dissipation ratio of the 3D Heisenberg SG 
based on off-equilibrium simulations [Hit ] . 

Now we turn to the diagonal-spin-overlap distribution 
function P(q dlag ) shown in Fig. 16(b). Although P(q d iag) 
exhibits a faint double-peak structure or the near flat- 
peak structure for smaller sizes of L < 12, it exhibits 
for larger sizes of L > 16 only a single peak located at 
qdiag = 0, which grows with increasing L, without any 
other appreciable peak structure. This is in contrast to 
the triple-peak structure observed in the chiral-overlap 
distribution function P(q x ) of Fig. 16(a), which is peaked 
at q x = and ±q GG . It is also in contrast to the double- 
peak structure observed in the spin-overlap distribution 
describing the ordered state of the mean-field Heisenberg 
SK model, which is peaked at q = ±\q EA (Hi- The ab- 
sence of any divergent peak at nonzero qdiag f° r larger 
L suggests that the model is in a SG disordered state, 
at least at a temperature T = 0.133. This conclusion 
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is consistent with our previous conclusion from the spin 
correlation-length ratio £,sg/L, the spin Binder ratio gsc 

(2) 

and the spin- glass order parameter q^Q, The appearance 
of a faint double-peak structure or a near flat-peak struc- 
ture for smaller sizes of L < 12 might be interpreted as 
a size-crossover from the small size SG pseudo-order to 
the large size SG disorder resulting from the expected 
spin-chirality coupling-decoupling behavior. 
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FIG. 16: (Color online) The overlap distribution function 
for the chirality (a), and for the spin (b), at a temperature 
T = 0.133. 

In Fig. 17, we show the temperature dependence of the 
non-self-averageness A parameters for the chirality (a) 
and for the spin (b) , respectively. Although the data are 
rather noisy due to large sample-to-sample fluctuations, 
the chiral Acg parameter of different L show a crossing 
and a prominent peak around the expected Tcg- This be- 
havior of Acg resembles the one observed in the 3D Ising 
SG with the Gaussian coupling (52j, the one of the 3D 
Heisenberg SG with the binary cou pling [l9j |. and the one 
of certain mean-field SG models 53,54]. At high temper- 
atures, both Acg and A$g tend to zero with increasing 
L, demonstrating that the system is self-averaging in this 
regime. Near T — Ted the peak hight of Aqg increases 
with increasing L, indicating that the system is non-self- 
averaging at Ted whereas below Tcg, the Acg still stays 
at nonzero value with increasing L, indicating that the 
CG ordered state is non-self-averaging. These findings, 
combined with the peculiar shape of P{q x ), suggest that 



the CG ordered phase accompanies an RSB with a non- 
self- averaging character. 

By contrast, the spin A$g parameter does not exhibit 
a peak at any temperature, but exhibits a crossing which 
occurs slightly above Tcg for the range of sizes stud- 
ied here. Below the crossing temperature, Asa tends to 
increase with L, suggesting that the XSG becomes non- 
self-averaging. Although the crossing of A$g is certainly 
a signature of a phase transition, it does not necessar- 
ily mean the occurrence of the standard SG transition 
characterized by a nonzero q SG or by the divergence of 
XSG- As mentioned in §3, a nonzero A$g persisting in 
the L — > oo limit simply means that the SG susceptibil- 
ity xsg is non-self-averaging. Below the CG transition 
temperature, one expects that the SG order parameter 
is still Gaussian distributed around zero with the width 
corresponding to a finite SG susceptibility xsGi while 
the width xsg/V~N exhibits samplc-to-sample fluctua- 
tions leading to the non-self- averaging (but finite) xsg- 
The latter is a natural consequence of the phase-space 
narrowing which should inevitably occur in the CG state 
exhibiting a one-step-like RSB. Of course, in the ther- 
modynamic limit, the width xsg/v^V vanishes yielding 
a 5-function located at g M „ = characteristic of the spin 
disordered state. Hence, the crossing of A$q curves and 
a nonzero- value of Asg remaining below Tcg are fully 
compatible with the absence of the standard SG long- 
range order below Tqq, which is consistent with our 
present observation of Tcg > Tsg ■ 

Fig. 18 exhibits the temperature dependence of the G 
parameters for the chirality (a) , and for the spin (b) . As 
can be seen from these figures, both the CG and SG G 
parameters exhibit a crossing near Tec, suggestive of a 
phase transition. As mentioned in §3, the G parameter 
can be written by the A parameter and the Binder ratio 
g as in Eqs. (20) and (21). As mentioned, a nonzero Gsg 
occurs once xsg becomes non-self-averaging, i.e., Gsg 

(2) 

0, even if there is no SG long-range order, i.e., qg G — 
and gsc = 0: See Eq.(21). Hence, the occurrence of 
a crossing in Gsg at Tec is entirely consistent with our 
observation of Tsg < Tec- We also note that the data of 
G are rather noisy with large error bars as compared with 
certain other quantities like the correlation-length ratio 
and the Binder ratio. Therefore, this quantity may not 
be well suited to an accurate estimate of the transition 
temperature T g . The same suggestion was also made by 
Ballesteros et al || and by Palassini et al [52j for the 
case of the 3D Ising SG. 



V. CRITICAL PROPERTIES OF THE 
CHIRAL-GLASS TRANSITION 

In this section, we study the critical properties of the 
CG transition on the basis of a finite-size scaling anal- 
ysis of our data of the CG susceptibility and the CG 
correlation-length ratio. 

From our analysis in the previous section, we fix the 
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FIG. 17: (Color online) The temperature and size dependence 
of the non-self-averageness A parameter for the chirality (a), 
and for the spin (b). An arrow indicates the location of the 
chiral-glass transition point. 



FIG. 18: (Color online) The temperature and size dependence 
of the G parameter for the chirality (a), and for the spin (b). 
An arrow indicates the location of the chiral-glass transition 
point. 



CG transition temperature to Tcg = 0.143 in this sec- 
tion. Our analysis in §4 already suggested the pres- 
ence of a significant correction-to-scaling term. Hence, 
we will try in this section to examine the effect of the 
correction-to-scaling, by setting the leading correction- 
to-scaling exponent to 9 = uj + i = 0.93 as was evalu- 
ated in §4. We then estimate the two independent criti- 
cal exponents characterizing the CG transition, i.e., the 
CG correlation-length exponent vcg and the CG critical- 
point-decay exponent rjcc- 

The standard finite-size scaling forms for the 
correlation-length ratio £,cg/L and for the CG suscep- 
tibility xcg are given by, 

^L=X((T-T CG )L 1 /^), (25) 

Xcg = L 2 -™°Y({T-T CG )L l ' v ° G ), (26) 

where X and Y are appropriate scaling functions. 

For the CG susceptibility xcg, we obtain a reasonably 
good scaling by the two-parameter fits with vcg — 1-3 ± 
0.2 and t)cg — 0.7 ± 0.2. The error bar quoted here and 
below is estimated by examining by eyes the quality of 
the fit with varying the fitting parameters. 



For the CG correlation-length ratio £,cg/L, our data 
shown in Fig. 9 have no common crossing point, indicating 
that the correction-to-scaling term is playing a significant 
role. Hence, we perform the scaling analysis of £cg/L 
with including the correction-to-scaling term, 

^ = X((T -Tc G )L 1/ucG )(l + aL-^), (27) 

where a is a numerical constant, and Tcg and u> are set 
Tcg = 0.143 and u + i = 0.93 as mentioned above. The 
resulting best scaling plot is shown in Fig. 19(a), to yield 
vcg — 1-4 ± 0.2. Note that by including the correction- 
to-scaling term we can obtain quite a good scaling, which 
has never been achieved unless we include the correction- 
to-scaling term in the analysis. This value of vcg — 
1.4 happens to be close to the one obtained from xcg 
without invoking the correction-to-scaling term. 

We also try a similar scaling analysis for xcc taking 
account of the correction-to-scaling term based on the 
form, 

Xcg = L 2 -^Y((T - T CG )L 1/UCG ){1 + aL-"). (28) 

The resulting best scaling plot is given in Fig. 19(b). The 
exponent estimates are vcg = 1.4 ±0.2 and rjcc = 0.6 ± 
0.2. 
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Recently, Campbell et al proposed an extended version 
of the standard finite-size scaling method, which might 
allow one to extend the scaling regime to a wider tem- 
perature range Q. In this method, one takes an appro- 
priate matching between the data in the critical regime 
and those in the higher temperature regime to extend 
the scaling regime. Campbell et al demonstrated that 
the method worked well for the 3D Ising SG Q. The 
relevant scaling forms are given by 

^=X((1-^)(LT) 1 ^), (29) 

Xcg = (LTf-^Y((l - Th-)(LTf'^% (30) 

f° r icc/L and xcg, respectively. 

Since this extended scaling method does not take care 
of the singular correction-to-scaling, it does not serve in 
itself to improve the quality of the scaling plot of £cg / L 
unless the singular correction-to-scaling term is invoked. 
If we apply the extended scaling form to xcg, we get a 
reasonably good scaling with vcg — 1.5 ±0.2 and rjcc = 
0.7 ±0.2, which are close to the values obtained based on 
the standard scaling form Eq.(26). 

It is also possible to apply this extended finite- 
size scaling both to £,cg/L and xcg with including 
the correction-to-scaling term. The appropriate scaling 
forms are given by 

^ = x((i - TM ){LT y/»c G){l + aL -^ (31) 

Xcg = {LTf-^oy ((1 _ Th.) {LT f /"co ){l + aL -u )t 

(32) 

The resulting best scaling plots are given in Fig. 20(a) for 
£cg/L, and in Fig. 20(b) for xcg- As can be seen from 
the figures, the quality of the fit is quite good, slightly 
better in a wider temperature region than the one ob- 
tained from the standard finite-size scaling. The expo- 
nent estimates are vcg — 1-5 ± 0.2 from ^cg/L, and 
v CG = 1.5 ± 0.2 and tj C g = 0.6 ± 0.2 from X cg- The 
value of vcg = 1-5 is slightly higher than the correspond- 
ing value vcg = 1-4 ± 0.2 obtained from the standard 
finite-size scaling of £cg/L, but they are fully compati- 
ble within the error bar. 

Combining the exponent estimates obtained from 
£,Cg/L and XCG, either by the standard analysis or by 
the extended one a la Campbell, we finally quote as our 
best estimates of the CG exponents, 

i/ OG = 1.4±0.2, tjcg = 0.6 ±0.2, (33) 

while the correction-to-scaling exponent takes a rather 
small value, u = 9 — - = 0.3 ± 0.1, suggesting that the 
correction-to-scaling effect is relatively large here. This 
value of lu is smaller than the corresponding value of the 
3D Ising SG, w ~ 1 gL H H. 
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FIG. 19: (Color online) Standard finite-size-scaling plots of 
the chiral-glass correlation- length ratio £cg/L (a), and of the 
chiral-glass susceptibility xcg (b), where the correction-to- 
scaling effect is taken into account. The chiral-glass transition 
temperature and the leading correction-to-scaling exponents 
are fixed to Tog = 0.143 and u + - = 0.93 as determined in 
§4. The best fit for £,cg/L is obtained with vcg = 1-4, while 
that for xcg is obtained with vcg = 1.4 and rjcc = 0.6. 



The estimated values of the CG critical exponents 
are compatible with the previous values obtained be- 
fore for the same model vcc — 1-2 and t\cg — 0-8 [ill ] 
and with those reported for the ±J 3D Heisenberg SG 
v CG = 1-2(2) and tj C g = 0.8(2) 0- By contrast, the 
obtained CG exponents differ significantly from the stan- 
dard exponent values of the 3D Ising SG, v ~ 2.5 ~ 2.7 
and 77 ~ -0.38 0.40 [1, Q. The result unambigu- 
ously indicates that the chiral-glass transition belongs to 
a universality class distinct from that of the standard 3D 
Ising SG, although the underlying Z2 symmetry is com- 
mon between the two. Possible long-range and/or many- 
body nature of the chirality-chirality interaction might 
be the cause of this difference. Further study is required 
to clarify the cause of this difference. 

Our scaling analysis in this section were based on the 
CG correlation-length ratio and the CG susceptibility. 
One may wonder if what happens if one uses the CG 
Binder ratio in the analysis. As is already evident from 
the form of gcc shown in Fig. 10(a), which exhibits a 
negative dip whose depth grows with the system size L, 
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FIG. 20: (Color online) Extended finite-size-scaling plots 
of the chiral-glass correlation- length ratio £cg/L (a), and 
of the chiral-glass susceptibility \cg (b) a la Campbell et 
al where the correction-to-scaling effect is taken into ac- 
count. The chiral-glass transition temperature and the lead- 
ing correction-to-scaling exponents are fixed to Tcg ~ 0.143 
and uj + i = 0.93 as determined in §4. The best fit for £cg / L 
is obtained with Vcg = 1-5, while that for xca is obtained 
with vcg = 1.5 and t\cg = 0.6. 



the finite-size scaling dose not work for gcG even with 
including the correction-to-scaling term. Such an exotic 
behavior of gec, e -<?-> the existence of a growing negative 
dip and the non-monotonic size dependence observed in 
certain temperature range above Tcg, is most probably 
reflecting the peculiarity of the CG ordered state itself, 
a possible one-step-like RSB feature, not just the sub- 
leading correction-to-scaling effect. If the peculiar behav- 
ior of gec arose reflecting the proximity of the nontrivial 
character of the CG ordered state, the finite-size-scaling 
analysis would not be applicable to gec m a straightfor- 
ward way, at least in the range of lattice sizes studied 
here. 

Finally, we wish to refer to the critical properties of the 
SG transition which is deduced to occur at T$g ~ 0.125. 
In §4, we already estimated the critical-point-decay expo- 
nent from the size dependence of the SG order parameter 

(2) <^ 

Isg as ^sg ~ —0.30. We also tried a fmite-size-scalmg 
analysis of both ^sg/L and xsd either with or with- 
out a correction-to-scaling term, just as we performed 



for the CG transition. However, it turns out that the 
finite-size scaling analysis for the spin-related quantities 
does not work well for any choice of T$g, v sg and tjsg, 
even if we adjust the assumed parameter values in a wide 
range. Possibly, the CG order that occurs preceding the 
SG order affects the scaling property of the spin-related 
quantities in a non-trivial way, and the data of still larger 
lattices and/or lower temperatures might be required to 
determine the critical properties of the SG transition. 



VI. SUMMARY 

In summary, we have studied equilibrium ordering 
properties of the three-dimensional isotropic Heisenberg 
spin glass by means of extensive Monte Carlo simulations. 
By calculating various physical quantities including the 
correlation-length ratio, the Binder ratio, the glass or- 
der parameter and the overlap distribution function up 
to the size as large as L = 32 and down to tempera- 
tures well below Tcg, we have given strong numerical 
evidence of successive CG and SG transitions occurring 
at T CG = 0.143±0.003 and at T SG < 0.125±g;g?|, respec- 
tively. The SG order sets in at a temperature at least 
about 10~15% below the CG order, hence, the occur- 
rence of the spin-chirality decoupling. On shorter length 
scale of L 5 12, the spin and the chirality often behave 
in a similar way, while, on longer length scale of L C 16, 
the chirality shows a stronger ordering tendency than 
the spin. The observation supports the view of the triv- 
ial spin-chirality coupling at shorter length scale crossing 
over to the spin-chirality decoupling at longer length scale 

MM- 

One may feel that the relative distance between Tcg 
and T$g is not so large, but, in fact, it is a sizable dif- 
ference, much larger than the one observed in other sys- 
tems exhibiting the spin-chirality decoupling, e.g., the 
2D regular frustrated XY model where the difference is 
known to be about 1% [H, US [HI . While the SG order 
in the 3D Heisenberg SG occurs at a nonzero temper- 
ature, as is consistent with the recent numerical works 
[13, |1 [H [H [H [H, H, HE HI, it should be stressed 
that whether T$g is zero or nonzero is irrelevant to the 
chirality scenario of Refs. l2ll [22J as long as the 

spin-chirality decoupling occurs, i.e., T$g < Tcg- 

We have observed a rather strong correction-to-scaling 
effect in our data, which, we have tried to control via the 
correction-to-scaling term with the correction-to-scaling 
exponent u> ~ 0.3. The analysis worked very well at least 
for the CG correlation length and the CG susceptibility. 
Hence, our conclusion of T$g < Tcg appears to be robust 
against the correction-to-scaling effect. 

We have also analyzed the critical properties associated 
with the CG transition. By mean of a finite-size scaling 
analysis with including the correction-to-scaling effect, 
we get an estimate of chiral-glass exponents vcg — 1.4 ± 
0.2 and ijcg — 0.6±0.2. The possibility of a simultaneous 
spin and chiral transition of the KT-type as suggested 
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in Refs.JHI, [36| is ruled out. We have shown that the 
behaviors of both the correlation-length ratio and the 
Binder ratio are entirely different from those of the 2D 
ferromagnetic XY model exhibiting the KT transition, 
and the KT scaling for these quantities does not work 
even with massive logarithmic corrections. 

The obtained values of the CG exponents are close to 
the values reported earlier in previous works, while they 
are entirely different from those of the 3D Ising SG. How- 
ever, these CG exponents are impressively close to the 
experimental values of SG exponents of canonical SGs 
like CuMn, AuFe, AgMn, etc, i.e., i/ ~ 1.3 ~ 1.4 and 
rj ~ 0.5 ~ 0.6. Indeed, this coincidence gives a strong 
support to the chirality scenario of experimental SG tran- 
sition of Ref.[l6|, [U, [2^] , since, in this scenario, 
the experimental SG exponents of weakly anisotropic 
Heisenberg-like SGs like canonical SGs are nothing but 
the CG exponents of the fully isotropic Heisenberg SG 
revealed via the random magnetic anisotropy. A very 
interesting consequence of the chirality scenario is that 
the chiral-glass transition, not the spin-glass transition, 
of the fully isotropic Heisenberg SG dictates the experi- 
mental SG transition. Experimentally, it remains highly 
interesting to directly estimate the set of chiral-glass ex- 
ponents by means of high-precision Hall measurements 
[59L l60l l6l| . It might also be worthwhile to re-examine 
the standard spin-glass exponents for various Heisenberg- 
like SG materials by controlling the magnitude of mag- 
netic anisotropy. 

Although various physical quantities have consistently 
suggested that the SG order occurs at a temperature 
lower than the CG transition temperature, a precise es- 
timate of the SG transition temperature and of the cor- 
responding SG exponents still remains to be a rather dif- 
ficult task, although we get t)sg ~ —0.30. In any case, 
the critical properties of the SG transition definitely dif- 
fer from those of the CG transition, since the associated 
r\ values are largely different. 

By measuring the Binder ratio, the overlap distri- 
bution function and the non-self-averageness parameter 
A, we have observed that the chiral-glass ordered state 
is non-self-averaging and exhibits a nontrivial phase- 
space structure (RSB). More precisely, we have observed 
a strong similarity to the systems exhibiting the so- 
called one-step RSB. We note that the one-step RSB 
feature was also observed in the same model in an off- 
equilibrium simulation probing the breaking pattern of 
the fluctuation-dissipation relation 15111 . According to 
the chirality scenario of Ref. [lfl [2(| l2ll . |22| , the proper- 
ties of the SG ordered state of real canonical SGs should 
be governed by the properties of the CG ordered state 
of the fully isotropic Heisenberg SG. If so, one-step-like 
RSB should eventually be an attribute of the SG ordered 
state of real canonical SGs. This is in sharp contrast to 
the long-standing common belief in the community, i.e., 
the SG ordered state of real canonical SGs exhibits either 
the hierarchical RSB (full RSB) or no RSB. 

After the submission of the manuscript, the authors 



learned that Fernadez et al also studied the same model 
by MC simulations up to the size L = 48, and suggested 
that the spin and the chirality might order simultane- 
ously [12]. We wish to give a few comments here: First, 
we have confirmed that the data of L < 32 reported in 
Ref. now agree with our present data within the er- 
ror bars. (This is somewhat in contrast to the data of 
Ref. [13| which deviate from our present data by 5 to 6 
of our a units.) Furthermore, the SG transition temper- 
ature reported in [62l | agrees with our present estimate 
Tsg = 0.125to o?2- The major difference then concerns 
with the difference in the estimate of the chiral-glass 
transition temperature Tqg- The reason of this discrep- 
ancy seems to be primarily originated from their L = 48 
£,cg/L data, which comes significantly smaller than the 
values expected from an extrapolation of the L < 32 data 
made in our present analysis. If the L — 48 chiral data of 
Ref. [1^] are to be trusted, it means a drastic changeover 
occurring in the chiral sector between the sizes L = 32 
and L = 48. The physical origin of this size crossover, 
if any, has yet to be identified. Meanwhile, since we do 
not have at the moment a plausible explanation of such a 
size-crossover, and since the equilibration of L = 48 chi- 
ral quantities is the hardest and the number of L = 48 
samples studied in Ref.[62| (164 samples) is significantly 
smaller than that of other sizes (984 samples), we feel 
that the L — 48 data of Ref. (62J should be cross-checked 
carefully by independent calculations. 

Overall, we believe that our present data give strong 
numerical support to the view that the spin-chirality de- 
coupling occurs in the 3D isotropic Heisenberg SG. Then, 
it also give support to the chirality scenario of experimen- 
tal SG transitions. 
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APPENDIX A: Simulation of the ferromagnetic 
2D XY model 

In this appendix, we report on the results of our simu- 
lation on the ferromagnetic XY (plane rotator) model on 
a 2D square lattice, a typical model exhibiting the KT 
transition. The Hamiltonian is given by 

n = -j^2 , (34) 

<ij> 

where J > is a ferromagnetic nearest-neighbor coupling 
and S{ is a two components classical unit vector at the 
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site i. The lattice is a Lx L square lattice (L ranging from 
8 to 512) with periodic boundary conditions. We perform 
an equilibrium MC simulation by using the single-spin- 
flip Metropolis method and the over-relaxation method. 
The over- relaxation sweeps are repeated M — L/8 times 
per every Metropolis sweep, which constitutes our unit 
MC step. 

The quantities we show here are the spin correlation- 
length ratio £/L and the spin Binder ratio g. The spin 
correlation-length is defined by 

, = i / MO) 2 ) " ^ 

* 2sin(fc m /2)\/ (m(fc m )2) W 

where k = (k m , 0) with fc m = 2tt/L, and (• • ■ ) denotes a 
thermal average, while 

1 N 

fi—x,y i—1 

is a ^-dependent magnetization, ./V being the total num- 
ber of the spins. The spin Binder ratio is defined by 
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FIG. 21: (Color online) The temperature and size dependence 
of the correlation-length ratio of the ferromagnetic 2D XY 
model. An arrow in the figure represents the location of the 
Kosterlitz-Thouless transition point. 

The temperature dependence of the correlation length 
ratio £/L is shown in Fig. 21. The KT transition tem- 
perature of this model was estimated rather precisely as 
T KT ~ 0.893 (in units of J) [&|. With increasing L, 
the £/L curves do not cross at a finite temperature, but 
tend to merge progressively at temperatures lower than 
Tkt, as can be seen from the figure. One can see that 
the observed behavior of £/L of the ferromagnetic 2D 
XY model is entirely different from the corresponding 
behavior of either the CG or SG correlation-length ra- 
tio, £,cg/L or (Isg/L, of the 3D Heisenberg SG shown in 



Fig. 9. The £,cg/L and £sg/L curves of the 3D Heisen- 
berg SG do not merge as in the £/L curves of the 2D XY 
model but intersect, the crossing points shifting to lower 
temperatures for large L 
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FIG. 22: (Color online) The temperature and size dependence 
of the Binder ratio of the ferromagnetic 2D XY model. An 
arrow in the figure represents the location of the Kosterlitz- 
Thouless transition point. 

Fig. 22 exhibits the temperature dependence of the 
Binder ratio g of the ferromagnetic 2D XY model. The 
g curves of different L now cross at a nonzero temper- 
ature. With increasing L, the crossing points approach 
Tkt from above. The data for larger L tend to give a 
"merging" behavior characteristic of the KT transition. 
This behavior is similar to the one reported by Loison for 
the same model Again, one sees that the observed 
behavior of g of the ferromagnetic 2D XY model is dif- 
ferent from the corresponding behavior of either the CG 
or SG Binder ratio, gcc or gsc, of the 3D Heisenberg 
SG shown in Fig. 10. 

Hence, from the comparison of the correlation-length 
ratio £ /L and the Binder ratio g of the ferromagnetic 2D 
XY model and of the 3D Heisenberg SG, one might also 
conclude that the transition of the 3D Heisenberg SG is 
not of the KT-type. 

APPENDIX B: Simulations of the ferromagnetic 
3D O(10) model 

In this appendix, we report on the results of our simu- 
lations on the ferromagnetic O(10) model on a 3D simple 
cubic lattice. The Hamiltonian is given by 

H = -J ^2 Si- Sj , (38) 

<ij> 

where J > is a ferromagnetic nearest-neighbor coupling 
and Si is a ten-components classical unit vector at the 
site i. In our simulation, we use the Mec algorithm of 
Ref.[65| combined with the over-relaxation method. The 
lattice is a L x L x L simple cubic lattice (L ranging from 
6 to 32) with periodic boundary conditions. The over- 
relaxation sweeps are repeated M = L/2 times per every 
Mec sweep, which constitutes our unit MC step. 
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As in Appendix A, the quantities we show here are 
the spin correlation-length ratio £/L and the spin Binder 
ratio g. The spin correlation- length ratio is defined by 



1 / (m(0) 2 ) 

2sin(fc m /2) y ( TO (A m ) 2 ) 



- 1, (39) 



where k = (k m , 0, 0) with k m = 2tt/L, and (• • • ) denotes 
a thermal average, while 



10 , N 



Mk) 2 = ^T\^J2S v eMik-n)\ 2 (40) 

fJ,— 1 2 — 1 

is a /c-dependent magnetization, N being the total num- 
ber of the spins. The spin Binder ratio is defined by 
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(m(0) 4 ) 
(m{6) 2 )' 



(41) 




FIG. 23: (Color online) The temperature and size dependence 
of the correlation-length ratio for the ferromagnetic 3D O(10) 
model. An arrow in the figure indicates the transition point. 
The inset exhibits the standard finite-size scaling plot without 
the correction term, where we put T c = 0.412 and v = 0.87. 

The temperature dependence of the correlation length 
ratio £/L is shown in Fig. 23. As can be seen from the 
figure, the £/L curves of various L show a clear crossing 
at an almost L-indepcndent temperature T c = 0.412 ± 
0.002 (in units of J), and splay out at lower temperatures. 

In Fig. 24, we show the temperature dependence of the 
Binder ratio g of the ferromagnetic 3D O(10) model. 



As can be seen from the figure, the g curves of differ- 
ent L also show a very clear crossing at an almost L- 
independent temperature T c = 0.412 ± 0.002, and splay 
out at lower temperatures. The behavior observed here 
for g is quite similar to the one observed for £/L in 
Fig. 23. In particular, in spite of its large number of 
order-parameter component of n = 10, no "merging" nor 
"marginal" behavior as suggested in Ref . [II] is observed. 
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FIG. 24: (Color online) The temperature and size dependence 
of the Binder ratio for the ferromagnetic 3D O(10) model. 
An arrow in the figure indicates the transition point. The 
inset exhibits the standard finite-size scaling plot without the 
correction term, where we put T c = 0.412 and v = 0.87. 



Very much similar "crossing" and "splaying out" behav- 
ior was observed also in the Binder ratio of the ferromag- 
netic 3D 0(6) model by Loison [50l | . 

This observation clearly demonstrates that the pecu- 
liar "non-crossing" and "negative dip" behavior as ob- 
served in the spin Binder ratio gsc of the 3D Heisenberg 
SG shown in Fig. 10(b) is not a trivial one originating 
from just the large number of order-parameter compo- 
nents (n = 9 in the case of the Heisenberg SG). As dis- 
cussed, the peculiar behavior observed in g$G of the 3D 
Heisenberg SG is likely to reflect an essential and peculiar 
feature of the ordered state of this model, most probably, 
the occurrence of a one-step- like RSB. 

Concomitantly, we also try to estimate the critical ex- 
ponent v from our data of £/L and g on the basis of the 
standard finite-size scaling analysis. Here, the correction 
seems to be negligible. Even without the correction term, 
we obtain a very good data collapse both for £/L and g, 
as shown in the insets of Figs. 23 and 24. We then get 
v= 0.87 ±0.03. 
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